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The successful operation of the relativistic heavy ion 
collider RHIC at BNL has stimulated an intensive study 
on the theoretically predicted quark gluon plasma (qgp) 
phase transition. It is expected that the appearance of 
phase transition will result in an anomalous behavior 
in the event-by-event fluctuations of transverse momen- 
tum iii. 

Since the multiplicity in a single event is finite, there 
are inevitably statistical fluctuations (sf) in the event-by- 
event analysis of transverse momentum. Various meth- 
ods Q ^ have been proposed to get rid of the statisti- 
cal fluctuations and measure the dynamical ones. Most 
of them are based on a subtraction procedure, i.e. first 
to estimate the variance of statistical fluctuations and 
then subtract it from that of the experimental data. The 
basic problem in this kind of methods is how to reliably 
estimate the statistical variance. 

Another approach proposed recently 7] is to eliminate 
the SF directly from the experimental data. This method 
can reproduce the dynamical variance exactly provided 
the SF are Poissonian, i.e. the particles are emitted in- 
dependently. 

The aim of this paper is to compare the different mea- 
sures proposed in the market and, in particular, to exam- 
ine the reliability of the estimation of statistical variance. 
The so called measure 8] will also be discussed and 
a new measure $p/ will be proposed, which can largely 
reduce the strong multiplicity dependence of itself. 

Let p{pt) be the dynamical transverse- momentum dis- 
tribution in single events and 



p{pt)dpt, (to= 1,2,...,M) (1) 



be the corresponding "coarse-grained" distribution 
The realization of p(pt) in experiment is the distribution 
of the total number n of particles in the pt region A, and 



the expression 

Qni = n,n/n (to = 1, 2, ... , M) (2) 

is an evaluation of , where n„j is the number of par- 
ticles falling into the mth bin. Thus the dynamical and 
experimentally-measured event-mean transverse momen- 
tum are, respectively. 



Ptc 



. M 
n = / PtP(Pt)dpt = V {Pt)mPm, (3) 
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where {pt)m is the pt value in the mth bin. The event- 
space moments of ptdyn ^^'^ Ptcxp ^^'^ 



CrP(Pt) = ((Ptexp)^) 




(5) 



(6) 



respectively. These two are evidently unequal due to 
the existence of SF. In particular, the experimental mo- 
ments C°^^{pt) depend crucially on the particle multi- 
plicity 11], which is mainly due to SF. Therefore, sub- 
traction procedures are commonly used to extract the 
dynamical fluctuation of pt . 

The basic idea of the subtraction procedure is that the 
variances <J^{pt) = C2{pt) — Ci(pt)^ of dynamical and 
statistical fluctuations are additive 
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(7) 
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where cr^^data: f^p-tdyn a^d cr^_j,tat are the experimental, 
dynamical and statistical variances, respectively. This 
additivity is true if the df and SF are independent, which 



2 



is a reasonable assumption. Therefore, the main problem 
in the subtraction procedure for extracting the dynami- 
cal variance is to estimate the statistical ones. Various 
methods have been proposed for this purpose. 

For example, in Ref. 4] the results from mixed events 
are considered as the baseline for the random distribu- 
tion and the difference in the fluctuation from a random 
distribution deflned as 



d — ^data ^baseline 



(8) 



is taken as a measure of the dynamical fluctuation. In 
Eq.(8) 



(pt) ' 



(9) 



Pt is the mean transverse momentum in a single event 
and (• • ■) denotes the average over event sample. 

Alternatively, in Ref. [3 the statistical variance of 
event mean pt , under the assumption of independent par- 
ticle production, is estimated as 
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ptstat 
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and, therefore, the dynamical variance is equal to 

2 2 2 



i^t VKR 



2 

ptstat 



Ptincl 



in) 



(10) 



(11) 



A somewhat different expression is used in Ref. ^Q] , where 
the experimental variance is weighted by event multiplic- 
ity A^,: 



Pt CER ES 



Ef=iN,{p-, 
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(12) 



where JV is the total number of events. 

A widely used measure for the non-statistical mean pt 
fluctuation is the proposed in Ref. jH 



(13) 



where z and Z are defined as z ; 



^ Pt — {ptj for each par- 
ticle and Z = y^ ''_i Zi = n{pt ~ {pt)) for each event, 
respectively |lCl||. The second term in the r.h.s. of 
Eq. (13) is the square root of the inclusive variance 
•^Ptincl = ((Pt - (Pt))^}- Assuming that the multiphc- 
ity fluctuation is uncorrelated with the pt fluctuation, we 
get from Eq.(13) 



V (n2)/(n)crp, 



data 



(T 



Ptincl- 



(14) 



This equation is evidently similar to Eq.(ll). 

All of the above measures have the same structure, be- 
ing based on a subtraction procedure, i.e. to subtract 
the variance of pt or a quantity related to it, that will be 
expected from a pure statistical system, from the same 
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FIG. 1: The experimental (open triangles) and dynamical 
(dashed line) variances of pt , together with the variances cal- 
culated from G-moments (open circles) and from Eqs. (11) 
(upward solid triangles) and (12) (downward solid triangles). 



quantity obtained in experiment. These measures will, 
of course, vanish for a pure statistical system, and a 
non-vanishing value of them will indicate the existence 
of dynamical effect. Therefore, the measures based on 
the subtraction procedure, as those listed above, will at 
least qualitatively measures the effect of df. 

In order to have an idea on how the various subtraction 
methods work quantitatively, let us carry on a Monte 
Carlo simulation using a toy model. In this model the 
distribution of pt is taken as 



p{Pt) = -^Pte 



-2pt/a 



(15) 



with a Gaussian distributed parameter a ((7^(a) — 0.24). 
In total 20 x 500, 000 events are generated for 20 average 
multiplicities. The resulting dynamical and experimental 
variances cr|-^dyn ^^"1 ""^texp ^'I'S plotted in Fig.l as dashed 
line and open triangles, respectively. For comparison, the 
CT^- vKR '''It CERES ^^'^ plotted in the same figure 
as upward and downward solid triangles, respectively. 
The differences between the experimental and dynam- 

are due to SF. It can be 



ical variances c^^exp 



cr. 



ptdyn 



seen from the figure that the cr^- vkr indeed reduces the 
SF, especially when the average multiplicity is not very 
low; while the ct^^ceres overshoots the influence of SF. 

On the other hand, the open circles shown in the figure 
reproduce the dynamical variances very well. They are 
obtained through a totally different approach 0, which 
is based on the elimination of the SF directly from the 
experimental data. 
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For the first order moment the ehmination of SF is provided the SF is Poissonian. The key point in eUminat- 



straightforward, 



M 



(16) 



ing SF in the second order moment is to expand the 2nd 
power m the definition Eq.(5) of C^'^^^^pt), 



Using the formulae 
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which holds for Poissonian SF, we then get Cj'^'*'"^ (pt) — G2{pt), where 



G2(Pt) = ( E(Pt)' 



\m— 1 



M 



E (Pt)m(Pt)r^ 



(19) 
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FIG. 2: Comparison of o^^^^^^^^ and crpj.i„;.i/{n) for systems 
with and without DF 



equality holds only for a pure statistical system without 
any dynamical fluctuation. When dynamical fluctuations 
are present the Cp- stat ^'^'^ '^ptinci/ ("■) ™ general, un 
equal. For an exhibition, the cr^ 



Ptstat 
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Fig.l, and cTp^ind/i^) plotted in Fig. 2 as sohd circles 
and solid triangles, respectively. They are approximately 
equal to each other only when the average multiplicity is 
not very low, e.g. when (n) > 12. 

In order to further clarify this point, we have made a 
pure statistical model with p„j — 1/M, cf. Eq.l. In this 
case crl-^dy^ = 0. The corresponding cr^^.t^ 

"'ptinci/i"') plotted in Fig. 2 as open circles and open 
triangles, respectively. The differences between them are 
small up to (n) = 5. This shows that the discrepancy 
between o-^_ vkr (o-p-tCEREs) and a^-^^^^ is mainly due to 



c^p-texp and 



Ptstat 



the inappropriate estimation at low multiplicity of dp- 
by o'pjinci/('^) for a system containing DF. 

Let us now turn to the discussion on the mea- 
sure The Monte Carlo results of are plotted 
in Fig. 3 as solid circles. It is strongly multiplicity depen- 
dent, and is, therefore, not a good measure for DF. On 
accounting of Eq.(14) we divide ^p^ by the square root 
of average multiplicity and define 



This guarantees the equality of G- variance ct'^^q = 
G2{Pt)-Gi{ptY to dynamical ones cr^- dyn C'2'^^"^(Pt)- 
Cf''"'(p-t)^ 

On the contrary, the SF cannot be exactly gotten rid of 
using the subtraction procedure. The main problem lies 
in the estimation of statistical variance by Eq.(lO). This 



$('■) = 
pt — 



cr. 



'^Ptdata 



Ptincl 



(20) 



cf. Eq.(ll). 

The Monte Carlo results of are shown in Fig. 3 
as upward solid triangles together with cr^t ceres ~ 

'^^tcEREs (downward triangles) and Op^^yi^ (dashed 
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FIG. 3: Variation of <E>pt and with (n), compared with 

'^pV CERES ^'^'^ '^Ptdyn 



line). It can be seen from the figure that has a much 



Ptdyn 



weaker multiphcity dependence and is closer to a 

than <I>pj . Therefore, is a better measure of df than 
<i>pj itself. 

In this paper, various measures for the dynamical 
transverse momentum fluctuation in event-by-event anal- 
ysis are compared with the dynamical fluctuation using 
a Monte Carlo model. 

It turns out that the ctJ ^ calculated from the G- 

Ptdyn' 



moments coincides with the dynamical variance 



showing that the G-moment method is effective in elimi- 
nating the statistical fluctuations coming from the finite 
number of particle in a single event and thus provides a 
good measure for the dynamical fluctuations. 

The cTp- VKR proposed by S. A. Voloshin, V. Koch and 
H. G. Ritter 5] and its revised version cr^^cERES M 
found to be good approximate measures for the dynam- 
ical fluctuation when the multiplicity is not very low. 
Therefore, their application to relativistic heavy ion ex- 
periments is justified. 

The divided by the square root of average multi- 
plicity (<I>p^ = ^pt/ \/ ('^)) is another approximate mea- 
sure for the dynamical fluctuation, which depends on 
multiplicity weaker than and is closer to the square 
root of dynamical variance CTp^dyn- 
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